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ABSTRACT 

We derive absolute dimensions of the early B-type detached eclipsing binary 
V453Cygni (B0.4IV + BO. 7 IV, P = 3.89 d), a member of the open cluster NGC 6871. 
From the analysis of new, high-resolution, spectroscopy and the UBV light curves 
of Cohen (1974) we find the masses to be 14.36 ± 0.20 M Q and 11.11 ± 0.13 M Q , the 
radii to be 8.55 ± 0.06 R© and 5.49 ± 0.06 Rq, and the effective temperatures to be 
26 600 ± 500 K and 25 500 ± 800 K for the primary and secondary stars, respectively. 
The surface gravity values of log g = 3.731 ± 0.012 and 4.005 ± 0.015 indicate that 
V453 Cyg is reaching the end of its main sequence lifetime. We have determined the 
apsidal motion period of the system to be 66.4 ± 1.8 yr using the technique of Lacy 
(1992) extended to include spectroscopic data as well as times of minimum light, giving 
a density concentration coefficient of logfc2 = —2.226 ± 0.024. Contaminating (third) 
light has been detected for the first time in the light curve of V453Cyg; previous 
analyses without this effect systematically underestimate the ratio of the radii of the 
two stars. The absolute dimensions of the system have been compared to the stellar 
evolution models of the Granada, Geneva, Padova and Cambridge groups. All model 
sets fit the data on V453 Cyg for solar helium and metal abundances and an age of 
10.0 ± 0.2 Myr. The Granada models also agree fully with the observed log/c2 once 
general relativistic effects have been accounted for. The Cambridge models with con- 
vective core overshooting fit V453Cyg better than those without. Given this success 
of the theoretical predictions, we briefly discuss which eclipsing binaries should be 
studied in order to further challenge the models. 

Key words: stars: binaries: eclipsing — open clusters — stars: fundamental pa- 
rameters — stars: binaries: spectroscopic — stars: early-type — stars: individual: 
V453Cygni 



1 INTRODUCTION 

Theoretical models of high-mass stars are very difficult to 
construct, due to a number of poorly understood astrophys- 
ical phenomena which become important at higher stellar 
masses. Convective core overshooting is an essential ingre- 
dient of stellar models and can have a large effect on the 
lifetimes and luminosities of high-mass stars. The degree of 
overshooting may depend on metallicity (Cordier et al. 2002, 
Palmieri etal. 2002) and stellar mass (Young etal. 2001). 
Stellar rotation is known to cause enhanced mixing in the 
interiors of stars which can mimic a small amount of convec- 
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tive overshooting and strongly affect the evolution of stars 
(Maeder & Meynet 2000). Mass loss is also an important 
effect (Maeder 1997) but is usually represented in theoret- 
ical models by a parameterisation such as that of Reimers 
(1975) (Woo etal. 2003). 

Until recently there was a discrepancy between masses 
derived from spectroscopic observations and inferred from 
evolutionary models (Herrero et al. 1992) . Recent advances 
in theoretical modelling, and the increased sophistication 
of observational techniques, appears to have removed this 
mass discrepancy (Hilditch, 2004), but confirmation of this 
requires detailed comparisons between models and as many 
physical properties of individual stars as possible. 

Detached eclipsing binaries (dEBs) are a vital source 
of observational data on the radii and masses of high-mass 
stars (Andersen 1991) but studies of such stars are hindered 
by the rarity of useful systems, the difficulty of measuring 
accurate velocities from spectra of rapidly-rotating multiple 



2 J. Southworth, P. F. L. Maxted and B. Smalley 



Table 1. Identifications, location, and combined photometric 
parameters for the V453 Cygni eclipsing system. 
References: (1) Cannon & Pickering (1923); (2) Argelander 
(1903); (3) H0g ct al. (1998); (4) Hoag etal.(1961); (5) Popper 
(1980); (6) Zakirov (1992); (7) Cohen (1969); (8) Reimann (1989). 





V453 Cygni 


Reference 


Henry Draper number 


HD 227696 


1 


Bonner Durchmusterung 


BD +35° 3964 


2 


Hoag number 


NGC 6871 13 


3 


"2000 


20 06 34.967 


4 


&200Q 


+35 44 26.28 


4 


Spectral type 


B 0.4 IV + B 0.7 IV 


5 


V 


8.285 


6 


B-V 


+0.179 


6 


U -B 


-0.61 


6 


V-R 


+0.254 


6 


P 


2.590 


7,8 



stars (Sana, Rauw & Gosset 2001), and the large amount 
of telescope time needed to obtain light curves which are 
complete over the whole orbital period. This has caused a 
shortage of data on eclipsing binaries with masses above 
10 Mq and sufficiently large orbital periods to be detached 
and therefore to have evolved as single stars. 

Whilst dEBs can provide accurate observed stellar 
masses, radii, effective temperatures and equatorial rota- 
tional velocities, the ages, chemical compositions and dis- 
tances of such systems are not, in general, directly observ- 
able. This means that when fitting observed quantities to 
theoretical predictions, there are several free parameters 
which can be adjusted to find the best fit to the observa- 
tional data. This makes it very difficult to assess the success 
of important but more subtle physical effects such as rota- 
tional mixing, mass loss, convective overshooting and even 
the type of turbulent convective theory. 

Eclipsing binaries in open clusters and associations can 
be analysed to derive accurate masses, radii, temperatures 
and rotational velocities. Their membership of a star cluster 
can also provide knowledge of the age, chemical composition, 
and distance of the dEB, allowing a more complete physical 
description of a single stellar system. This reduces the num- 
ber of free parameters adjustable when fitting theoretical 
models to observations, which allows a discriminate test of 
the representation of different physical effects and indeed the 
success of different sets of stellar evolutionary models. For 
example, from the analysis of two main sequence dEBs which 
are members of the open cluster h Persei (NGC 869) , South- 
worth, Maxted & Smalley (2004, hereafter Paper I) were able 
to provide a first estimate of the bulk metallicity of one of 
the most important and well-studied open clusters in the 
Northern Hemisphere. 



1.1 V453 Cygni in NGC 6871 

V453 Cygni is a high-mass dEB with an orbital period 
of 3.89 days. Its membership of the young open cluster 
NGC 6871 means that its age and distance can be found 



independently. The primary component of V453 Cyg is ap- 
proaching the terminal age main sequence (TAMS) and its 
large radius causes the eclipses to be total, allowing a very 
accurate determination of the radii of both stars. Table 
contains identifications and some photometric properties of 
the system. 

The eclipsing nature of V453 Cyg was discovered by 
Wachmann (1939) and an early spectroscopic orbit was cal- 
culated by Pearce (1941). A period study by Cohen (1971) 
provided a determination of the orbital longitude of perias- 
tron, lu, inconsistent with that derived by Pearce. In a period 
study by Wachmann (1973) this was correctly interpreted 
as apsidal motion. Wachmann derived an apsidal period of 
U = 72 years using measurements of the time differences 
between several groups of adjacent primary and secondary 
eclipses. A more recent period study, including parabolic and 
periodic terms, was undertaken by Rafert (1982). 

Excellent photoelectric UBV light curves were observed 
by Wachmann (1974) and analysed using the Russell-Merrill 
method (Russell & Merrill 1952) which involves the proce- 
dure of rectification. His work contains a plot of the light 
curves adjusted for the effects of orbital eccentricity and 
apsidal motion using parameters updated from that of his 
previous work, but the data themselves have so far been 
unobtainable. It is possible that they are in an unlabelled 
file in the IAU Variable Star Archives (Breger 1988), but no 
record of them exists at Hamburg Observatory, where the 
light curves were observed (A. Reiners, private communica- 
tion). 

Cohen (1974) published complete photoelectric UBV 
light curves which contain fewer datapoints and more obser- 
vational scatter than those of Wachmann (1974). He anal- 
ysed these using the Russell-Merrill method but stated that 
his observations were not definitive. They have since been 
analysed by Cester et al. (1978) using the light curve analysis 
code wink (Wood 1971). This is the only previous photo- 
metric study to use modern techniques. 

A recent investigation using photoelectric UBVRI light 
curves has been published by Zakirov (1992). He analysed 
his light curves using the "direct machine method of Lavrov 
(1993)", which is based on rectification. The results of the 
four photometric analyses of V453 Cyg are substantially in 
agreement about the basic photometric parameters of the 
system. 

Recent spectroscopic orbits have been published by 
Popper & Hill (1991), Simon & Sturm (1994) and 
Burkholder, Massey & Morrell (1997, hereafter BMM97). 
These results are collected in Table Simon & Sturm used 
seven spectra to demonstrate their spectral disentangling 
algorithm, which decomposes observed composite spectra 
into the separate spectra of two stars. The total secondary 
eclipse of V453 Cyg allowed them to directly compare their 
disentangled primary spectrum with a spectrum observed 
during secondary eclipse. Disentangling can be used to de- 
termine accurate orbital semiamplitudes (Hynes & Maxted 
1998, Harries et al. 2003) but it is not clear if there is a 
robust method by which to estimate uncertainties in the de- 
rived quantities. BMM97 derived a good spectroscopic orbit 
from 25 high signal-to-noise spectra and compared the dEB 
to models to investigate the discrepancy at higher masses 
between models and observations. The rotational velocities 
of the components of V453 Cyg were determined by Olson 
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Table 2. Published spectroscopic orbits of V453Cygni. BMM97 originally fitted their data with a circular orbit. We have refitted their 
radial velocities with an eccentric orbit to increase the accuracy of our determination of the apsidal motion (see section|3]for details). A 
colon after a number indicates that it is uncertain. Quantities without quoted errors or a colon were not determined by that investigation. 
When quantities are given separately for each star we have quoted a weighted mean of the two values. Symbols have their usual meanings. 
Times are written as (HJD - 2 400 000). 

* The reference time, To, refers to a time of periastron passage, not a time of minimum light. 





Pearce 


Abt, Levy and 


Popper and 


Simon and 


BMM97 


BMM97 




(1941) 


Gandet (1972) 


Hill (1991) 


Sturm (1994) 




(our solution) 


P (days) 


3.87972 


3.8890 


3.8898128 


3.88982309 


3.8898128 


3.889825 


T (HJD) 


30231.0843 ±0.0543* 


40495.027* 


39340.099 


36811.7296 


48141.82 ±0.01* 


48500.64 ±0.66* 


K A (kms" 1 ) 


181. 8± 1.13 


152: 


171 ±1.5 


171.7±2.9 


173.2 ±1.3 


173.7 ±1.4 


K B (kms" 1 ) 


237.4 ±2.78 




222 ±2.5 


223.1 ±2.9 


213.6 ±3.0 


212.4 ±3.4 


e 


0.07 ±0.007 


0.05: 


0.0 


0.0 


0.0 


0.011 ±0.015 


lu (degrees) 


175.2 ±5.06 


99: 








88.6 ±6.0 


V 1 (kms" 1 ) 


-15.0±0.94 


-22.7: 


-14 


-7: 


-17.6 ±1.0 


-18.0±1.6 



(1984) to be 107 ±9 kms" 1 and 97 ±20 kms" 1 for the pri- 
mary and secondary stars respectively. A preliminary single- 
lined spectroscopic orbit was also given by Abt, Levy & 
Gandet (1972). 

An abundance analysis of V453 Cyg was undertaken by 
Daflon et al. (2001) using both LTE and non-LTE calcula- 
tions. The results suggest that V453 Cyg has a slightly sub- 
solar metallicity. These authors derived an effective temper- 
ature of 29 100 K using the Q parameter based on U BV mag- 
nitudes (Johnson 1957), and a surface gravity of log g = 4.45 
(cms -2 ) from profile fitting of the H7 4340 A spectral line. 
Both values are larger than expected and inconsistent with 
previous analyses. Their surface gravity value is in fact some- 
what larger than appropriate for the ZAMS, and is incon- 
sistent with significant main sequence evolution. 

I. 2 NGC 6871 

The open cluster NGC 6871 is a concentration of bright OB 
stars which forms the nucleus of the CygOB3 association 
(Garmany & Stencel 1992). This makes it an important ob- 
ject for the study of the evolution of high-mass stars. The 
cluster itself has been studied photometrically several times 
but its sparse nature means determination of its physical 
parameters is difficult. 

UBV photometry of the 30 brightest stars was pub- 
lished by Hoag etal. (1961). Crawford, Barnes & Warren 
(1974) observed 11 stars in the Stromgren uvby system and 
24 stars in the Crawford /3 system, finding significantly vari- 
able reddening and a distance modulus of 11.5 mag. This 
uvby/3 photometry was extended to 40 stars by Reimann 
(1989), who found reddening E),- y with a mean value of 
0.348 mag and an intracluster variation of about 0.1 mag. 
His derived distance modulus of 11.94 ±0.08 mag and age of 
12 Myr are both greater than previous literature values. 

Massey, Johnson & DeGioia- Eastwood (1995) con- 
ducted extensive UBV CCD photometry of 1955 stars in 
the area of CygOB3. Their values of distance modulus, 

II. 65±0.07, and of reddening, E B -v = 0.46±0.03 mag with 
individual values between 0.04 and 1.11 mag, agree well with 
previous determinations. They find an age of 2 to 5 Myr for 
stars with spectral types earlier than B0 but give evidence 
for a significant spread of stellar ages in the cluster. Whilst 



the highest-mass unevolved cluster members have main se- 
quence lifetimes of 4 to 5 Myr, NGC 6871 contains evolved 
15 Mq stars despite their main sequence lifetimes being of 
the order of 11 Myr. 



2 OBSERVATIONS 

Spectroscopic observations were carried out in 2002 October 
on the 2.5 m Isaac Newton Telescope (INT) on La Palma. 
The 500 mm camera of the Intermediate Dispersion Spec- 
trograph (IDS) was used with a holographic 2400 1mm" 1 
grating and EEV 4k x 2k CCD. From measurements of the 
full width half maximum (FWHM) of arc lines taken for 
wavelength calibration we estimate that the resolution is 
approximately 0.2 A. The spectral windows chosen for ob- 
servation were 4450-4715 A (31 spectra) and 4230-4500 A 
(12 spectra). Additional spectra were observed around H/3 
(4861 A) to provide an additional temperature indicator for 
spectral analysis. The signal-to-noise ratio per pixel of the 
observed spectra is between 100 and 450. 

Data reduction was undertaken using optimal extrac- 
tion as implemented in the software tools PAMELA and 
molly 1 (Marsh 1989). 



3 PERIOD DETERMINATION AND APSIDAL 
MOTION 

A period study by Wachmann (1973, 1974) indicated fast 
apsidal motion with a period of U — 71 yr from eighteen 
times of minimum light covering almost three thousand or- 
bital cycles. Photographic minima exist dating back to the 
year 1902 but they are not of sufficient quality to improve the 
apsidal motion analysis (Ashbrook, unpublished, but tabu- 
lated in Cohen 1971). Zakirov (1991) states that his obser- 
vations are not consistent with Wachmann's apsidal motion 
period. 

Times of minima for inclusion in the apsidal motion 
analysis were taken from Cohen (1971), Wachmann (1973), 

1 PAMELA and MOLLY were written by Dr. Tom Marsh and can 
be found at http://www.warwick.ac.Uk/staff/T.R.Marsh 
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Table 5. Apdisal motion parameters for V453 Cyg. To illustrate the final result we have also included the best-fitting parameters using 
only photometric data, and the results of Wachmann (1974). Times are written as (HJD — 2 400 000). 





Wachmann 


This paper 


This paper 




(1974) 


(photometric only) 


(final results) 


Anomalistic period P (days) 




3.890426 ±0.000073 


3.890450 ±0.000017 


Reference minimum time To 


36811.7296 


39340.1011 ±0.0022 


39340.0998 ±0.0019 


Orbital inclination i (degrees) 




88.0 (fixed) 


88.0 (fixed) 


Orbital eccentricity e 


0.02 


0.019 ±0.002 


0.022 ±0.002 


Periastron longitude at To, uiq (degrees) 


309.2 


313.2 ±6.7 


309.7 ±3.1 


Apsidal motion rate u> (degrees P s ~~ 1 ) 


0.0539 


0.0556 ±0.0068 


0.0579 ±0.0016 


Sidereal period P s (days) 


3.88982309 


3.889824 ±0.000082 


3.889825 ±0.000018 


Apsidal motion period U (years) 


71 


68.9 ±8.5 


66.4 ±1.8 
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Figure 1. Representation of the best-fitting apsidal motion parameters. The upper panel shows the observed times of primary (circles) 
and secondary (triangles) minima, minus the expected times given by a linear ephemeris, compared to the best-fitting curves of primary 
(dashed) and secondary (dotted) minima. The open circle represents the rejected time of minimum of Biro etal. (1998). The lower panel 
shows the spcctroscopically determined longitudes of periastron, ui, and the change of ui over orbital cycle. Errors have only been shown 
if they arc larger than the corresponding symbol. 



R. Diethelm 2 , and Biro etal. (1998). The method of Lacy 
(1992) was adopted to solve the apsidal motion equations 
without the use of approximations. Here the Levenberg- 
Marquardt non-linear least-squares fitting routine (Press 
et al. 1992) is used to relate the derived apsidal motion pa- 
rameters to the observed times of minima. The authors' im- 
plementation of this method (apsmot) uses subroutines very 



2 Eclipsing Binaries Minima Database at 
http : //www . oa . uj . edu . pl/ktt /index . html 



generously supplied by D. Holmgren (see Holmgren & Wolf 
1996). 

It was immediately clear that the more recent times of 
minima were not in full agreement with each other or with 
the times of minima used by Wachmann (1973). In such 
cases, independent information is needed to decide which 
published observations are reliable and which arc discrepant. 
For this reason we added to our APSMOT code the ability to 
include spectroscopic determinations of eccentricity, e, and 
longitude of periastron, w, in the overall fit. 

Values of e and lu were taken from the spectroscopic 
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Table 3. Times of minimum light of V453 Cyg taken from the 
literature. The O — C values refer to the difference between the 
observed and calculated values. 
* Rejected from the fit due to a large O — C value. 
References: (1) Wachmann (1973) photographic, (2) Wachmann 
(1973) photoelectric, (3) Cohen (1971) photoelectric, (4) R. Di- 
ethelm (see text) photoelectric, (5) Biro etal. (1998) CCD. 



Cycle Minimum time Adopted O — C Ref. 
number (HJD-2 400 000) error 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
2 
2 
2 
3 
3 
3 
3 
2 
2 
4 
5 



Table 4. Spectroscopic data used in the apsidal motion analysis. 
References: (1) Pearce (1941), (2) BMM97 (our solution). 



Cycle Eccentricity O-C tv O-C Ref. 

number (e) (degrees) (w) 

-2342.0 0.070 ±0.007 0.048 175.2 ±5.1 1.0 1 

2355.0 0.011 ±0.015 0.011 88.6 ±6.0 2.6 2 



studies by Pearce (1941) and BMM97. The radial velocity 
observations of BMM97 were originally fitted with a circular 
spectroscopic orbit so we have reanalysed the velocities (Ta- 
ble |5| to determine e and u>, and assigned a cycle number 
corresponding to the approximate midpoint of the observa- 
tions. We have fitted an eccentric orbit to the data using 
SBOP 3 in single-lined Lehman-Filhes mode. There is a large 
correlation between ui and the ephemeris reference time, To, 
which makes both values somewhat uncertain. The solution 
of the secondary velocities did not converge without fixing 
the value of To to that of the primary star's spectroscopic 
orbit. 

The spectroscopic data of Pearce (1941) and BMM97 
allowed us to identify the time of minimum of Bi'ro etal. 
(1998) as being in disagreement with the other photometric 

3 Spectroscopic Binary Orbit Program written by Dr. Paul B. 
Etzel (http : //mintaka. sdsu. edu/f aculty/etzel/). 



Table 6. Equivalent widths of helium lines in the spectra of 
V453Cyg. These are true equivalent widths per individual star, 
after corrections for dilution due to the spectra being composite. 



Species 


Wavelength 


Equivalent widths (A) 




of line (A) 


primary star 


secondary star 


He I 


4387.93 


0.661 


0.890 


He I 


4437.55 


0.089 




He I 


4471.50 


0.992 


1.18 


Hen 


4685.70 


0.139 


0.057 



data. This datapoint has been rejected from the apsidal mo- 
tion solution, and the other data were assigned appropriate 
uncertainties. These data are given in Tables H and H along 
with the assigned uncertainties and the observed minus cal- 
culated (O — C) values. The final solution is plotted against 
the data in Fig. Q and given in Table where it is com- 
pared to the solution of Wachmann (1974) and to a solution 
without the inclusion of spectroscopic data. 

The time of minimum of Bi'ro et al. (1998) has been con- 
firmed by reanalysis and another unpublished time of min- 
imum (I. Biro, private communication). If they are correct, 
they may indicate the existence on another effect on the 
times of minima, for example the light-time effect. Further 
data are needed to investigate this possibility. 



4 SPECTRAL SYNTHESIS 

The observed spectra were fitted to synthetic spectra, by 
the method of least squares, to derive the effective tempera- 
tures of the components of V453 Cyg. Synthetic spectra were 
calculated using uclsyn (Smith 1992, Smalley etal. 2001), 
Kurucz (1993) ATLAS9 model atmospheres and absorption 
lines from the Kurucz & Bell (1995) linelist. The profiles 
for the 4387.93 A and 4471.50 A Hei lines were calculated 
using profiles from the work of Barnard etal. (1969) and 
Shamey (1969), with log gf values from the critical compi- 
lation of Wiese etal. (1966). The spectra were rotationally 
broadened as necessary and an instrumental broadening of 
FWHM = 0.2 A was applied to match the resolution of the 
observations. 

The primary star was analysed using a spectrum ob- 
tained during a total secondary eclipse. For the secondary 
star we used the spectra at quadrature to measure the 
lines and corrected them for dilution effects. The equiva- 
lent widths of the helium lines are given in Table inland were 
used to obtain values of effective temperature (T e ff) by en- 
suring ionisation balance between He i and He ii, for assumed 
values of surface gravity and microturbulence velocity (£t). 

Using the surface gravities found from the spectroscopic 
and photometric analyses in this work, we find the effective 
temperatures T e ft = 26 600 ± 500 K for the primary star 
and T c ff = 25 500 ± 800 K for the secondary star. These 
parameters imply that the atmospheres of these stars are 
helium-rich by about 0.25 dex compared to solar. Further 
support for these T e s and logg values is given by the H7 
4340 A profiles; H7 profiles with higher logg (as found by 
Daflon et al. 2001) are too broad to fit the observations. 

Using the uvby/3 photometry from Hauck & Mermilliod 



-2790.0 


28487.531 


0.01 


0.0026 


-2789.5 


28489.435 


0.01 


0.0008 


-2608.0 


29195.476 


0.01 


-0.0028 


-2607.5 


29197.371 


0.01 


-0.0018 


-1501.0 


33501.508 


0.01 


-0.0020 


-1500.5 


33503.414 


0.01 


-0.0008 


-1482.0 


33575.411 


0.01 


-0.0054 


-1481.5 


33577.328 


0.01 


0.0061 


-1390.0 


33933.270 


0.01 


-0.0084 


-1389.5 


33935.195 


0.01 


0.0074 


-65.0 


39087.266 


0.005 


0.0039 


-64.5 


39089.242 


0.005 


0.0028 


-7.5 


39310.9552 


0.005 


-0.0052 


-7.0 


39312.8702 


0.005 


-0.0010 


-5.5 


39318.7347 


0.005 


-0.0054 


-5.0 


39320.6497 


0.005 


-0.0011 


12.0 


39386.7764 


0.005 


-0.0010 


178.0 


40032.492 


0.005 


0.0068 


178.5 


40034.470 


0.005 


-0.0015 


1684.0 


45890.5660 


0.005 


0.0016 


2801.0 


50235.4843* 


0.005 


-0.0426 
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Table 7. Radial velocities and O — C values (in kms" 1 ) for 
V453 Cyg calculated using TODCOR. Weights are given in column 
"Wt" and were derived from the amount of light collected in that 
observation and were used in the SBOP analysis. 



Table 8. Parameters of the spectroscopic orbit derived fom 
V453 Cyg using TODCOR only on narrow lines. The systemic veloc- 
ities were derived using TODCOR and synthetic template spectra. 



HJD - 
2 400 000 


Primary 
velocity 


O-C 


Secondary 
velocity 


O-C 


Wt 


52562.5002 


-188.9 


-0.4 


201.3 


-9.3 


1.0 


52564.3199 


142.9 


-4.0 


-218.9 


4.2 


1.0 


52566.2961 


-181.4 


1.8 


200.7 


-3.0 


0.7 


52566.3112 


-191.0 


-6.7 


204.5 


-0.7 


0.7 


52566.4356 


-191.5 


-1.9 


190.6 


-21.4 


1.2 


52568.3471 


150.8 


-5.2 


-227.8 


7.0 


1.0 


52560.3187 


137.6 


3.5 


-206.5 


0.1 


1.0 


52560.3210 


136.0 


1.5 


-211.7 


-4.6 


1.0 


52560.3233 


133.0 


-1.8 


-217.9 


-10.5 


1.0 


52560.3257 


133.8 


-1.3 


-218.2 


-10.3 


1.0 


52560.3280 


131.1 


-4.3 


-211.7 


-3.4 


1.1 


52560.4091 


142.5 


-2.4 


-219.0 


1.1 


0.8 


52560.4114 


146.4 


1.2 


-221.5 


-0.7 


0.8 


52560.4137 


149.3 


4.0 


-225.9 


-4.8 


0.8 


52560.4161 


144.7 


-0.9 


-212.8 


8.6 


0.8 


52560.4184 


151.3 


5.4 


-225.9 


-4.2 


0.8 


52562.4949 


-184.9 


3.4 


219.0 


8.6 


1.2 


52564.3111 


147.4 


1.3 


-217.4 


1.7 


1.1 


52566.3018 


-189.0 


-5.4 


216.5 


12.3 


1.1 


52566.3165 


-183.5 


1.1 


216.5 


10.9 


1.0 


52568.3558 


160.0 


3.7 


-227.5 


7.7 


1.3 


52570.3545 


-184.4 


5.4 


216.5 


1.1 


1.3 



(1998) and the uvbybeta and TEFFLOGG codes of Moon 
(1985), we have obtained de-reddened photometry and fits 
to the grids of Moon & Dworetsky (1984). Values of T cS = 
26 710 ± 800 K and logg = 3.78 ± 0.07 were obtained, in 
excellent agreement with the parameters of the primary star, 
which dominates the system light, obtained in this work. 

Daflon et al. (2001) adopted T cft = 29 100 K, log g = 
4.45 (cms -2 ) and £t = 12 kms -1 in their detailed analysis 
of V453 Cyg. The above ionisation balance analysis gives 
a T cff = 29 200 K for their logg and While this is in 
agreement with the T e g they adopted, their value of logg 
is not supported by our absolute stellar dimensions, our H7 
profile fitting and the uvbyP photometry, so we prefer cooler 
effective temperatures for the components of V453 Cyg. 



5 SPECTROSCOPIC ORBITS 

Radial velocities of the two stars were derived from the ob- 
served spectra using the two-dimensional cross-correlation 
algorithm TODCOR (Zucker & Mazeh 1994). In this method 
the observed spectra are cross-correlated against two tem- 
plate spectra simultaneously. Whilst in Paper I we used syn- 
thetic spectra from uclsyn as templates, here we can use 
spectra obtained around the midpoint of secondary eclipse. 
As the eclipses are total, this contains only light from the 
primary star (and a negligible amount of contaminating light 
- see section [fjj. This template was used for both stars due 
to the similarity of their spectral characteristics, and allows 
us to avoid possible systematic errors due to the use of theo- 
retical spectra. Best-fitting synthetic spectra were also gen- 



Primary Secondary 


Orbital period P (days) 


3.889825 (fixed) 


Reference time To (HJD) 


39340.6765 (fixed) 


Eccentricity e 


0.022 (fixed) 


Periastron longitude w (°) 


140.1 (fixed) 


Semiamplitude K (kms -1 ) 


173.7 ±0.8 224.6 ±2.0 


Systemic velocity (kms -1 ) 


-13.1 ±0.3 -16.2 ±1.8 


Mass ratio q 


0.773 ± 0.008 


a sin i ( Rq ) 


30.59 ± 0.17 


M sin 3 i (M Q ) 


14.35 ±0.20 11.10 ±0.13 



200 r 
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Figure 2. Spectroscopic orbit for V453 Cyg from an SBOP fit to 
radial velocities from TODCOR. 



erated and an independent solution was performed using 
TODCOR to check our observed-template solution and deter- 
mine the systemic velocities of the stars. 

TODCOR was unreliable when the velocity separation of 
the stars was significantly lower than the combined rota- 
tional velocities of the stars. For this reason the best spectra 
were selected by eye, leaving sixteen covering the wavelength 
region of 4450-4715 A and six spectra covering 4230-4500 A. 
The hydrogen and helium lines at 4340 A, 4471 A and 4686 A 
were masked to avoid significant errors due to the blending 
of broad spectral lines (Andersen 1975). The resulting ra- 
dial velocities are given in Table and were fitted using 
SBOP in single-lined Lehman-Filhes mode. The orbital pe- 
riod, ephemeris time of reference, eccentricity and ui were 
fixed at values derived from the apsidal motion analysis. 
The results of this analysis are given in Table GO and the 
final spectroscopic orbit is plotted in Fig. [5] 

Our velocity semiamplitudes are slightly larger, al- 
though in general consistent with, those found in the recent 
spectroscopic analyses of Popper & Hill (1991), Simon & 
Sturm (1994) and BMM97. This effect is probably because 
our radial velocities have been derived using only metal lines, 
whereas previous studies have relied mainly on helium lines. 
The effect of neglecting orbital eccentricity, however, is neg- 
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ligible, as can be seen from the two orbital solutions of the 
BMM97 radial velocities in Tabled 



6 LIGHT CURVE ANALYSIS 

We have analysed the U BV light curves taken from the work 
of Cohen (1974). As discussed in section FTTI these observa- 
tions are not definitive, but for this totally eclipsing system 
they are able to provide accurate values of the individual 
stellar radii. We have used the simple and efficient NDE 
eclipsing binary model as implemented in the light curve 
analysis program ebop 4 (Nelson & Davis 1972, Popper & 
Etzel 1981). In this model the stellar shapes are approxi- 
mated by biaxial ellipsoids. The calculated oblatenesses of 
the best-fitting model for V453 Cyg are within the limits of 
reliability for the ebop code (Popper & Etzel 1981). 

The light curves were phased with the sidereal period. 
Filter-specific linear limb darkening coefficients (LDCs) were 
taken from van Hamme (1993), gravity darkening exponents 
pi were fixed at 1.0 (Claret 1998) and the mass ratio was 
fixed at the spectroscopic value. After an initial solution was 
obtained, datapoints which showed residuals of more than 
3a were rejected. These observations were clearly discrepant 
and their omission has not affected the derived parameter 
values but has lowered their uncertainties. 

Initial investigation suggested that there is a small 
amount of third light, L3, but acceptable solutions can be 
found without this effect. However, solutions with L3 7^ fit 
slightly better than solutions with L3 = for all three light 
curves, so third light has been included in all final solutions. 
If Z/3 is neglected, a ratio of the radii lower by about 0.04 is 
required to reproduce the observed eclipse depths. The ef- 
fect on the derived stellar radii is an increase in Ra by about 
1% and a decrease in Rb by about 4%. These adjustments 
would bring our photometric solution into agreement with 
previous light curve analyses, which have all neglected third 
light and therefore may be systematically wrong. 

Table [5] shows several theoretical determinations of lin- 
ear LDCs for stars having similar effective temperatures and 
suface gravities of the components of V453 Cyg. We have 
evaluated the effect of a change in LDCs, on the parameters 
of the photometric solution. Table |U] suggests that there is 
a variation of about 0.05 between different investigations of 
LDCs, so we perturbed the van Hamme (1993) values by 
this amount and refitted the light curves. The resulting er- 
rors have been added to the quoted uncertainties in Table HUl 
but are significant only for the surface brightness ratios. We 
have also determined the upper values of the LDCs for which 
light curve fits are not notably worse than our best fits, as- 
suming that the limb darkening of both stars are adequately 
represented by the same LDC. This shows that the recent 
trend towards larger theoretical LDCs is producing values 
which are close to the observational maxima for our B and 
V light curve solutions. 

The best-fitting light curves are compared to the ob- 
servations in Fig. [3] The residuals of the fit are also shown, 
and some minor systematic trends are noticeable. Whilst the 

4 Eclipsing Binary Orbit Program written by Dr. Paul B. Etzel 
(http://mintaka.sdsu.edu/faculty/etzel/). 



ebop light curve model is adequate to fit the current photo- 
metric data, definitive light curves may require a more so- 
phisticated treatment such as that contained in the Wilson- 
Devinney code (Wilson & Devinney 1971, Wilson 1993), 
which has a better representation of limb darkening and the 
reflection effect. V453 Cyg is a good system for the deter- 
mination of observational LDCs due to the long totality of 
its primary eclipse. The Cohen (1974) light curves are not 
of sufficient quality to determine LDCs here; definitive light 
curves will be required. 

6.1 Error analysis 

We have modified the ebop code to use the Levenberg- 
Marquardt minimisation algorithm (mrqmin; Press et al. 
1992) to find the best fit to the data. Whilst mrqmin al- 
lows calculation of the formal errors of the adjusted light 
curve parameters, it is known that these uncertainties are 
very optimistic when some parameters are significantly cor- 
related. Correlations are generally small for systems which 
exhibit total eclipses, except for systems with third light. 
Orbital inclination and L3 can be strongly anticorrelated 
as both have a significant dependence on the depth of the 
eclipses. Robust estimation of uncertainties must include an 
assessment of parameter correlations for the physical char- 
acteristics of the system under investigation. 

We have used the method of bootstrapping to evaluate 
the uncertainties and correlations of the parameters derived 
from the light curve analysis. After the best fit was deter- 
mined for each light curve, a synthetic light curve was eval- 
uated at the phases of observation of the real light curve. 
We added observational noise of the same magnitude as the 
real light curve and refitted the synthetic light curve with 
ebop. This process was undertaken ten thousand times for 
each observed light curve. 

It is important to understand what information these 
bootstrapping simulations actually provide. Once a best fit 
is found, the distributions of the ten thousand evaluations 
of various parameters give us the parameter uncertainties 
and their correlations, based on the best fit, the phases of 
observation, and the observational scatter of the real light 
curves. This is a valid method of analysis if the best fits are 
close to the true characteristics of the dEB. The reality of 
this assumption can be assessed using independent solutions 
of different light curves, for example the U, B and V obser- 
vations here. The bootstrapping analysis then serves as an 
indication of the validity of uncertainties estimated from the 
interagreement of different light curves. 

Sample plots of the distributions of different parameter 
values are shown in Fig. [I] It is notable that the ratio of 
the radii, k, and the ratio of the surface brightnesses, J, are 
not correlated due to this system exhibiting total eclipses. 
However, k and L3 show a very strong correlation and illus- 
trate why L3 has not been included in previous light curve 
analyses. This effect is because, for a given value of k, a 
well-defined value of L3 is required to fit the well-determined 
eclipse depths. 

The best-fitting photometric parameters of V453 Cyg, 
their la uncertainties, and the final adopted parameters are 
given in Table HOI The adopted parameters were determined 
using weighted means of the values determined from the 
individual light curves, and their uncertainties are similar 



8 J. Southworth, P. F. L. Maxted and B. Smalley 



Table 9. Theoretical limb darkening coefficients (LDCs) for stars similar to the components of V453 Cyg, compared to constraints from 
analysis of the light curves. A subscripted A or B denotes the LDC of the primary star or secondary star respectively. Filter designations 
are given in brackets. 



Reference 


u A (U) 


u B (U) 


u A (B) 


u B (B) 


u A (V) 


u B (V) 


Klinglesmith & Sobieski (1970) 


0.340 


0.321 


0.317 


0.302 


0.268 


0.249 


Wade & Rucinski (1985) 


0.327 


0.318 


0.292 


0.279 


0.274 


0.248 


vanHamme (1993) 


0.324 


0.296 


0.318 


0.287 


0.260 


0.247 


Di'az-Cordoves, Claret & Gimenez (1995) 


0.374 


0.357 


0.376 


0.357 


0.334 


0.316 


Claret (1998) 


0.370 


0.360 


0.365 


0.358 


0.320 


0.309 


Claret (2000) 


0.423 


0.384 


0.420 


0.381 


0.374 


0.334 


Largest LDCs which fit the light curves well 


0.5 


0.4 


0.35 



Table 10. Results of the light curve analysis for V453 Cygni. The adopted values are the weighted means of the values determined from 
the individual light curves. 







U 


B 


V 


Adopted 


Total number of datapoints 




538 


540 


540 


1618 


Number used in solution 




531 


532 


534 


1597 


Linear limb darkening coefficient 


UA 


0.324 


0.318 


0.260 




Linear limb darkening coefficient 


"B 


0.296 


0.287 


0.247 




Primary radius (a) 


VA 


0.2788 ±0.0021 


0.2800 ±0.0014 


0.2793 ±0.0014 


0.2795 ±0.0009 


Secondary radius (a) 


'B 


0.1781 ±0.0039 


0.1785 ±0.0030 


0.1811 ±0.0029 


0.1794 ±0.0018 


Ratio of the radii 


k 


0.648 ±0.016 


0.637 ±0.013 


0.649 ±0.012 


0.644 ±0.008 


Orbital inclination (degrees) 


% 


89.9 ±1.3 


88.2 ±1.2 


89.0 ±1.1 


89.0 ±0.7 


Surface brightness ratio 


J 


0.938 ±0.014 


0.953 ±0.010 


0.948 ±0.015 




Light ratio 


Lb/ La 


0.381 ±0.022 


0.375 ±0.016 


0.384 ±0.015 




Third light (fraction of total light) 


L 3 


0.079 ±0.025 


0.068 ±0.024 


0.089 ±0.022 





to but larger than the standard deviations of the individual 
values. 



7 ABSOLUTE DIMENSIONS AND 

COMPARISON WITH STELLAR MODELS 

The derived physical parameters for the component stars of 
V453 Cyg have been collected in Table 1111 the masses and 
radii of the two stars have been measured to accuracies bet- 
ter than 1.4%. The radius of the primary star is extremely 
well determined because it depends mainly on the duration 
of totality during secondary eclipse. The rotational velocities 
of the stars are slightly uncertain: using our spectroscopic 
data we have been unable to derive values more accurate 
than those of Olson (1984). The primary star rotates syn- 
chronously with the orbital velocity but the secondary ro- 
tates somewhat faster (although with a large uncertainty). 

The timescale for orbital circularisation is about 1.5 Gyr 
(Zahn 1977; Hilditch 2001) which is consistent with the pres- 
ence of a small but well-defined eccentricity. The timescale 
for synchronisation is approximately 3.0 Myr (Zahn 1977; 
Hilditch 2001) and the rotational velocity of the secondary 
star suggests that V453 Cyg is younger than this. As we 
have measured the age of V453 Cyg to be 10.0 ± 0.2 Myr by 
comparison with the predictions from theoretical models, we 
expect that this slight discrepancy is due to the secondary 
rotational velocity being quite uncertain. 



7.1 Stellar model fits 

The absolute parameters of the components of V453 Cyg 
have been compared to the predictions of stellar models 
from four different groups:- (1) Claret (1995), hereafter re- 
ferred to as the Granada95 models, (2) Bressan et al. (1993), 
hereafter referred to as the Padova93 models (the more 
recent Padova models of Girardi etal. (2000) only extend 
to stellar masses of 7M©), (3) Schaller etal. (1992), here- 
after referred to as the Geneva92 models, and (4) Pols et al. 
(1998) with modifications 5 , hereafter referred to as the Cam- 
bridge2000 models. For each set of models we have interpo- 
lated over age using cubic spline functions and plotted the 
resulting predictions in the mass-radius and logarithmic ef- 
fective temperature-surface gravity diagrams. Comparisons 
with the properties of V453 Cyg were performed simultane- 
ously in both diagrams and the two stars were assumed to 
have the same age and chemical composition (as expected for 
close binary stars). The Granada95, Padova93 and Geneva92 
models all include a moderate amount of convective core 
overshooting (although with different formalisations) . Hap- 
pily, the Cambridge2000 models are available both with and 
without a moderate amount of overshooting, allowing us to 
test whether the inclusion of this effect provides a better fit 
to the observational data. 



For details of these modifications please refer to the website of 
O. R. Pols: http://www.astro.uu.nl/~pols/ 
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Figure 3. Observed phased light curves of V453 Cyg with the best-fitting EBOP model light curves. The lower three curves show the 
residuals of the EBOP fits. For clarity the B and V residuals are offset by +0.15 and +0.3 magnitudes, respectively. 



Panels (a) and (b) of Fig. show the parameters of 
V453 Cyg compared to the predictions of the Granada95 
models. A good fit is obtained for an age of 9.9 Myr and 
a chemical composition of (Z,Y) = (0.02,0.28) (i.e., normal 
helium abundance). Attempts to fit the stars with a higher 
or lower helium abundance (Claret 1995) or metal abun- 
dances of Z = 0.01 (Claret & Gimenez 1995) or Z = 0.03 
(Claret 1997) were unsuccessful. 

The predictions of the Padova93 and the Geneva92 
models are compared to the parameters of V453 Cyg in pan- 
els (c) and (d) of Fig. |^1 The two sets of model predictions 
are plotted for Z — 0.02 and the same three ages, so are 
directly comparable apart from a slight difference in the 
assumed helium abundance (Y = 0.28 for the Padova93 
models and Y = 0.30 for the Geneva92 models). It is no- 
table that the two sets of models agree very well not only 
with each other but with the comparable Granada95 models 



discussed above. Whilst both the Padova93 and Geneva92 
models fit the components of V453 Cyg best for an age of 
9.8 Myr, a marginally better fit is provided by the slightly 
higher predicted T e $ values of the Padova93 models. At- 
tempts were also made to fit the components of V453 Cyg 
using metal abundances higher [Z = 0.05 for the Padova93 
models (Fagotto etal. 1994a); Z = 0.04 for the Geneva92 
models (Schaerer etal. 1993b)] and lower [Z = 0.008 for 
the Padova93 models (Fagotto etal. 1994b); Z = 0.008 for 
the Geneva92 models (Schaerer etal. 1993a)] without any 
success. 

The Cambridge2000 model set differs from the other 
model sets considered here in that it is available with and 
without a moderate amount of convective core overshoot- 
ing, but does not include any mass loss; this should be 
unimportant for these stars. Panels (e) and (f) of Fig. 
show the parameters of V453 Cyg compared to the predic- 
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Figure 5. Comparison between stellar models and the absolute dimensions of V453 Cyg in the mass-radius and the T e ff-log g diagrams. 
Isochrones have been plotted to represent the stellar models, with circles showing their points of evaluation. Broken lines have been 
plotted by interpolating over mass using cubic splines. Panels (a) and (b) show the Granada stellar models for (X,Y) = (0.02,0.28) and 
ages of 9.7, 9.9, 10.1 and 10.3 Myr (radii increase and T c ffS decrease as age increases). Panels (c) and (d) show the Padova models for 
(X,Y) = (0.02,0.28) (dotted lines) and Geneva models for (X,Y) = (0.02,0.30) (dashed lines) for ages of 9.4, 9.8 and 10.2 Myr. Panels 
(e) and (f) show the Cambridge models for (X,Y) = (0.02,0.28) with overshooting (dashed lines) and without overshooting (dotted lines, 
for ages of 9.8 and 10.2 Myr (with overshooting) or 9.4 and 9.8 Myr (with no overshooting). 



Eclipsing binaries in open clusters. II. V453 Cyg in NGC6871 11 




coot 

0.50 



0.70 0.50 0.60 0.70 0.50 0.60 0.70 



0.70 0.50 



0.70 0.50 



0.60 0.70 



Table 11. Absolute dimensions of the detached eclipsing binary 
V453Cygni in the open cluster NGC6871. 

* Calculated using the combined magnitude and flux ratio in the 
V filter, the assumed cluster distance modulus and reddening, 
and the canonical reddening law Ay = 3.1Eb—v- 

' Calculated using the effective temperature— bolometric correc- 
tion calibration of Flower (1996) 

* Taken from Olson (1984). 
V e q and Kjynch refer to the 
tional velocities respectively. 



Cluster age logr (years) 
Cluster distance modulus 



V453 Cyg A V453 Cyg B 

6.3 to 6.7 
11.65 ± 0.07 



Orbital period (days) 


3.889825 ± 0.000017 


Mass ratio q 


0.773 ± 0.008 


Mass (Mq) 


14.36 ±0.20 


11.11 ±0.13 


Radius ( R0) 


8.551 ±0.055 


5.489 ±0.063 


log g ( cm s — 2 ) 


3.731 ±0.012 


4.005 ±0.015 


Effective temperature (K) 


26 600 ±500 


25 500 ± 800 


My* (mag) 


-4.44 ±0.38 


-3.39 ±0.39 


Luminosity t (log L/Lq) 


4.69 ±0.21 


4.24 ±0.28 


V cq t (kms" 1 ) 


107 ±9 


97 ±20 



Kynch (kms" 1 ) 111.3 ±0.7 71.4 ±0.8 

Systemic velocity (kms -1 ) -13.1±0.3 -16.2±1.8 



Apsidal motion period (yr) 
log k 2 



66.4 ± 1.8 
-2.226 ± 0.024 



4Q0 E 1 ' L 1 3 E 1 1 3 E 1 : 1 

0.50 0.60 0.70 0.50 0.60 0.70 0.50 0.60 0.70 




35 56 97 38 39 35 35 57 3 3 59 35 96 37 36 89 
inclination inclination inclination 



Figure 4. Sample distributions of the best-fitting parameters 
evaluated during the bootstrapping analysis. The units and pa- 
rameter symbols are as in Table Hol Each distribution between two 
parameters is shown for the U (left), B (middle) and V (right) 
light curves. 

tions of the Cambridge2000 models. Both overshooting and 
standard-mixing isochrones are plotted for an age of 9.8 Myr 
for comparison. The overshooting models are also plotted 
for the best-fitting age of 10.2 Myr and the standard-mixing 
models are plotted for their best-fitting age of 9.4 Myr. The 
overshooting models are notably more successful than the 
standard-mixing models, which predict T c g values which 
are slightly too low for the stars of V453 Cyg. As with the 
Padova93 and Geneva92 models, we were unable to perform 
fits for different helium abundances as such models have not 
been published. 

The above comparisons demonstrate that a good agree- 
ment has been reached between different sets of theoretical 
evolutionary models for stars similar to the components of 
V453 Cyg. All sets of models were successful in fitting the 



observations for an age of 10.0 ± 0.2 Myr and solar metal 
and helium abundances. We also attempted to fit models to 
the absolute dimensions of V453 Cyg derived with zero third 
light. This changes the radii to 8.649 and 5.250 R© and the 
surface gravities to logy = 3.723 and 4.045 (cms" 2 ), with 
other quantities, and the uncertainties, unaffected. Using the 
Granada95 models we were able to achieve a fit in the mass- 
radius plane for low metal abundance (Z — 0.01), high he- 
lium abundance (Y — 0.36) and an age of 8.2 Myr. However, 
the T c ff values were predicted to be 2000 K greater than ob- 
served, and a combination of low metal abundance and high 
helium abundance does not agree with the predictions of 
Galactic chemical evolution theory (see e.g., Binney & Mer- 
ryfield 1998). We were unable to find a simultaneous fit in 
the mass-radius and T c g-\ogg diagrams for the Geneva92, 
Padova93 or Cambridge2000 models. 

BMM97 successfully fitted the Geneva92 theoretical 
models to the observed masses and luminosities of the com- 
ponents of V453 Cyg. As they did not compare stellar radii, 
they were not subject to errors from the assumption of no 
third light. The luminosities of V453 Cyg are also much 
more uncertain than the radii derived here, so fitting in the 
log T c ff — Mboi plane allows a wider range of predictions to 
fit the observed data. 

7.2 Comparison between the observed apsidal 
motion constant and theoretical predictions 

The internal structure constant log &2 is a measure of the 
density concentration of components of a dEB which shows 
apsidal motion. As discussed in Claret & Gimenez (1993), 
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the observed density concentration coefficient can be calcu- 
lated from the apsidal motion period using the equation 

1 P 

C21 + C 2 2 U 

where P and U are the orbital period and apsidal period 
and the constants C2; are weights which depend on the stel- 
lar mass, radius and rotational velocity (i=l refers to the 
primary star and i=2 refers to the secondary) . The individ- 
ual theoretical density concentration coefficients, ft2i, must 
be combined using the equation 



C21 + c 2 2 

to find the weighted average coefficient which is directly 
comparable the observed value. 

We must remember that the observed log k-z contains 
not only the Newtonian contribution to apsidal motion, 
which is due to the stars being distorted extended shapes, 
but a small amount of general relativistic apsidal motion. 
This last contribution can be straightforwardly calculated 
using 

LU = 5.45 X 10 t; 

1 - e 2 V P ) 

(Gimenez 1985) where uj is the apsidal motion rate in degrees 
per cycle, e is the orbital eccentricity, P is the anomalistic 
period (days) and Mi are the stellar masses (Mq). 

The tabulated results for the Granada95 theoretical 
models include predicted internal structure constants which 
can be compared with our observations. The general rela- 
tivistic contribution to the apsidal motion of V453 Cyg is 
approximately 6% of the total apsidal motion rate, and its 
subtraction gives a Newtonian density concentration coeffi- 
cient of 

log fc 2 obs = -2.254 ± 0.024 

(neglecting any uncertainty in calculation of the general rela- 
tivistic contribution). The theoretical weighted average coef- 
ficient was calculated using individual fe; values from inter- 
polation, in the table of model predictions, to the observed 
masses and best-fitting age of V453 Cyg. This value is dom- 
inated by the primary star, as it is more distorted due to its 
significantly larger radius, and is 

log fc 2 tho ° = -2.255 

The agreement with observations is excellent. This agree- 
ment is particularly important for assessing the assumed 
amount of convective core overshooting, on which theoreti- 
cal values of log k2 have a significant dependence (Claret & 
Gimenez 1991). 

8 MEMBERSHIP OF THE OPEN CLUSTER 
NGC 6871 

V453 Cyg is traditionally considered to be a member of the 
NGC 6871 open cluster, and appears on the cluster main se- 
quence in all photometric diagrams. Further proof of mem- 
bership comes from its systemic velocity, — 13.2 ± 0.3 km s _1 
using a weighted mean of systemic velocities calculated for 
each star. This agrees well with the value of —15 ± 6kms _1 
quoted by Hron (1987) and the radial velocity of the cluster 



member NGC 6871 13 which was measured to be —14.6 ± 
2.7 km s -1 using the same instrumental setup as we used for 
V453 Cyg. However, the radial velocity of NGC 6871 is given 
as -7.7i3.2kms- 1 by Rastorguev etal. (1999), which dif- 
fers from our value for V453 Cyg by 1.6a. The proper motion 
of V453 Cyg is consistent with cluster membership (Perry- 
man etal. 1997). 

Massey etal. (1995) give an age of 2 to 5Myr for 
the members of NGC 6871 with the earliest spectral types, 
but their photometric diagrams contain somewhat evolved 
15 M Q stars which are also claimed to be cluster members. 
This suggests that the stars in NGC 6871 have either a 
spread in ages or were created by two distinct bursts of star 
formation. We cannot currently distinguish between the two 
possibilities; the age of 10.0 ± 0.2 Myr derived for V453 Cyg 
using theoretical models is consistent with both evolutionary 
scenarios. 



9 DISCUSSION 

We have derived the absolute dimensions of the components 
of the high-mass detached eclipsing binary V453 Cygni, a 
member of the open cluster NGC 6871. 

Effective temperatures were found using the helium ion- 
isation balance derived from high-resolution spectra, which 
also suggest an enhanced photospheric helium abundance 
relative to solar. Radial velocities were derived from the 
spectra using only the weak spectral lines and the TODCOR 
two-dimensional cross-correlation algorithm. 

The apsidal motion rate of the system has been deter- 
mined using an extended version of the photometric method 
of Lacy (1992), which includes times of minimum light and 
spectroscopic determinations of eccentricity and ui. The ap- 
sidal period is well constrained, and allow the derivation of 
eccentricity and w to a greater accuracy than possible with 
the light curves and radial velocity curves. 

We have reanalysed the UBV light curves of Cohen 
(1974) in order to determine the radii of the components 
of the dEB. The best-fitting parameters include a small 
amount of third light, which was previously undetected. Ro- 
bust parameter uncertainties were derived using bootstrap- 
ping, allowing us to quantify and illustrate the effect of cor- 
relations between different photometric parameters. The ra- 
tio of the radii and the amount of third light are strongly 
correlated, due to their dependence on the depths of the 
eclipses; previous photometric studies which did not include 
third light are systematically biased towards values of the 
stellar radii which are 1% higher and 5% lower for primary 
and secondary respectively. 

The accurate absolute dimensions presented here allow 
V453Cyg to be added to the list of dEBs with the best- 
determined values of mass, radius and effective temperature 
(Andersen 1991). However, our analysis would clearly be 
much improved with better observational data. The inclu- 
sion of only a few new times of minima would greatly in- 
crease the accuracy of the results of the apsidal motion anal- 
ysis, and more accurate rotational velocities would allow a 
more accurate derivation of the internal structure constant 
logfc2- A definitive spectroscopic orbit will require obser- 
vations with a higher signal-to-noise than those presented 
here, and should give masses determined to accuracies of 
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better than 1%. Definitive light curves of the system would 
allow determination of the limb darkening coefficients for 
both stars, providing an important test of model atmosphere 
codes. 

The absolute masses, radii and effective temperatures of 
the components of V453 Cyg have been compared to several 
stellar models in the mass-radius and log T e ff-log g planes, 
assuming the same age for both stars. Not only is there im- 
pressive agreement betwen different theoretical models, all 
model sets are able to fit the observational data for a solar 
helium and metal abundance. Moreover, the Granada mod- 
els (Claret 1995) provide a perfect match to the observed ap- 
sidal motion rate once the relativistic contribution has been 
subtracted from the overall effect. Stellar models have for 
a long time appeared to predict that the central condensa- 
tions of stars are lower than that found using observations 
of apsidal motion (e.g., Clausen, Gimenez & Scarfe 1986, 
Barembaum & Etzel 1995). This apparent discrepancy has 
been reduced by the discovery that the internal structure 
constants change significantly through a star's evolution. 
The current generation of theoretical models, incorporating 
OPAL opacity data (Rogers & Iglesias 1992), are in good 
agreement with observations. It is noticeable that some ob- 
servers have not removed the general relativistic effect from 
their observed log &2 values before comparison with theory; 
in many cases this will have a negligible effect but for the 
stars of V453 Cyg it causes about 6% of the observed ap- 
sidal motion, changing log by an amount similar to its 
uncertainty. 

The normal helium abundance implied by stellar model 
fits also conflicts with the slight overabundance noted in our 
spectral synthesis analysis. We note that the photospheric 
helium abundance is not directly comparable to the initial 
internal helium abundance used in model calculations. 

Fits to the Cambridge stellar models support the inclu- 
sion of a moderate amount of overshooting in most stellar 
evolutionary models. Whilst models without overshooting 
were able to fit the masses and radii of the stars, the pre- 
dicted effective temperatures are slightly lower than that 
determined from the helium ionisation balance. 

The stellar models were extremely successful in fitting 
the absolute dimensions and effective temperatures of a 
high- mass slightly-evolved dEB, with component masses and 
radii differing by ten and twenty-five times their combined 
uncertainties, respectively. For observational stellar astro- 
physicists, this fact implies that we must either observe sys- 
tems so thoroughly that their masses and radii are known 
to accuracies of 0.5% and the effective temperatures to 2%, 
or target particular types of stars to critique the success 
of one set of stellar models compared to another. Such tar- 
gets include low-mass, high-mass, pulsating, and Population 
II stars, as well as eclipsing systems found in Local Group 
galaxies. Eclipsing binaries in open clusters can satisfy this 
requirement if the cluster they belong to is well-studied or 
otherwise interesting (Paper I), and further observations are 
being undertaken towards this goal. 
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